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ABSTRACT 

We address a long-standing discrepancy between non-LTE analyses of the prominent C n AA4267 and 
6578/82 A multiplets in early- type stars. A comprehensive non-LTE model atom of Cn is constructed 
based on critically selected atomic data. This model atom is used for an abundance study of six 
apparently slow-rotating main-sequence and giant early B-type stars. High- resolution and high-S/N 
spectra allow us to derive highly consistent abundances not only from the classical features but also 
from up to 18 further C n lines in the visual - including two so far unreported emission features 
equally well reproduced in non-LTE. These results require the stellar atmospheric parameters to be 
determined with care. A homogeneous (slightly) sub-solar present-day carbon abundance from young 
stars in the solar vicinity (in associations and in the field) of logC/H+12 = 8.29 ± 0.03 is indicated. 
Subject headings: line: formation - radiative transfer - stars: abundances, early type 



1. INTRODUCTION 

One of the most abundant metals in the universe is car- 
bon, the central building block of all organic chemistry. 
Abundances derived from luminous early B- and late O- 
type stars provide important constraints on stellar and 
galactochemical evolution. In extragalactic applications 
(e.g. analyses of dwarf stars in the Magellanic Clouds, 
see e.g. Hunter et al. 2005; Korn et al. 2005) one desires 
to study the strongest features in the metal line spectra, 
as these are the only measurable at low S/N and/or high 
projected rotational velocities. In the case of ionized car- 
bon these are the prominent lines of the two multiplets 
Cn AA4267.02/4267.27A and 6578.03/6582.85 A. These 
lines are unfortunately highly sensitive to non-LTE ef- 
fects, as well as to the choice of stellar atmospheric pa- 
rameters. So far, all studies from the literature failed to 
derive consistent abundances from these lines. 

The problem was addressed by Lennon (1983), Eber 
& Butler (1988) and Sigut (1996: S96) using Cn model 
atoms of increasing complexity. Spectrum synthesis 
based on the latter two non-LTE model atoms obtained 
better - but apparently still not good - agreement with 
observation when compared to earlier work or LTE re- 
sults. Non-LTE abundance analyses in Galactic OB stars 
e.g. by Gies & Lambert (1992) andKilian (1992) employ- 
ing the Eber & Butler (1988) model derived a metal de- 
ficiency (including C) in these young stars with respect 
to the Sun, in accordance with studies of H n regions. 

The present work aims at providing a solution to the 
classical non-LTE problem of carbon abundance deter- 
minations in OB stars. A reliable C II model atom is de- 
veloped and first applications on high-quality spectra are 
presented. Besides great care in the selection of atomic 
data, special emphasis is also given to an accurate atmo- 
spheric parameter determination, both in order to mini- 
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mize systematic uncertainties. The measurable C II spec- 
trum in the visual is investigated, the two prominent line 
multiplets as well as numerous weaker lines. 

2. MODEL CALCULATIONS 

A hybrid approach is used for the non-LTE line forma- 
tion computations. These are based on line-blanketed 
LTE model atmospheres calculated with ATLAS9 (Ku- 
rucz 1993). Synthetic spectra are computed with re- 
cent versions of DETAIL and SURFACE (Giddings 1981; 
Butler & Giddings 1985), solving the restricted non-LTE 
problem - see Przybilla et al. (2001) for details. 

Hydrogen and He I /n populations are computed using 
recent model atoms by Przybilla & Butler (2004) and 
Przybilla (2005), respectively. We have compared our 
synthetic H and He I /n lines with predictions from Lanz 
& Hubeny (2003) in the effective temperature range of 
27500 K < T e g < 32500 K for dwarf as well as giant stars. 
Overall good agreement is found. Exceptions are the Hei 
singlet lines, which are predicted (significantly) stronger 
in our approach, in accordance with observation. 

The C II model atom considers LS-terms up to prin- 
cipal quantum number n = 10 and angular momentum 
£ = 9 explicitly in non-LTE, with all fine-structure sub- 
levels combined into one. The doublet and quartet spin- 
systems are treated simultaneously. Level energies are 
adopted from Moore (1993), S96 and Quinet (1998). 

Oscillator strengths from three sources are considered: 
fine-structure data from ab-initio computations using the 
multiconfiguration Hartree-Fock method in the Breit- 
Pauli approximation (Froese Fischer & Tachiev 2004: 
FFT04), data from an application of the Breit-Pauli R- 
matrix method (Nahar 2002a) and results obtained in 
the Opacity Project (OP) from a close-coupling method 
in the LS-approximation (Yan, Taylor & Seaton 1987). 
Our primary source of (//-values is FFT04, followed by 
OP and Nahar for the remaining transitions. Note that 
data for several important transitions from Nahar dis- 
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agree with those from FFT04 and OP, which show consis- 
tency among each other. Intercombinations are neglected 
because of the very small ^/-values and high densties. 

Photoionization cross-sections are taken from the OP 
(Yan & Seaton 1987) where available. For the remainder, 
data from Nahar (2002b) are adopted. The choice is em- 
pirically motivated, giving preference to the OP data in 
order to reproduce observation over the entire parameter 
range simultaneously from all indicators. The two data 
sets show differences in the resonance structures. 

Effective collision strengths for electron impact exci- 
tation among the lowest 16 LS-states are adopted from 
R- matrix computations of Wilson, Bell & Hudson (2005). 
Collisional excitation for transitions without reliable 
data are treated using the Van Regemorter (1962) ap- 
proximation in the optically allowed case and via the 
semi-empirical Allen (1973) formula otherwise. Colli- 
sion strengths varying between 0.01 (An > 4) to 100 
(An = 0) are employed, as suggested by evaluation of 
the detailed data of Wilson et al. (2005). 

Collisional ionization rates are evaluated according to 
the Seaton (1962) approximation, with threshold pho- 
toionization cross-sections from OP and Nahar (2002b), 
allowing for an empirical correction of one order of mag- 
nitude higher for the 6f 2 F° and 6g 2 G terms - corre- 
sponding to the upper levels of the Cn AA6151 and 
6462 A transitions. For completeness, a C in model atom 
is also accounted for in the computations, but its details 
are of no further importance to the present study. 

Voigt profiles are adopted in the formal solution using 
SURFACE. Wavelengths and (//-values are taken from 
Wiese et al. (1996). Radiative damping parameters are 
calculated from OP lifetimes and coefficients for colli- 
sional broadening by electron impact are adopted from 
Griem (1974) for the Cn A4267A doublet, while the ap- 
proximation of Cowley (1971) is used for the other lines. 

3. ANALYSIS 

Six apparently slow-rotating stars are considered for 
the model atom calibration and first applications. The 
observations consist of high-resolution, high-S/N spectra 
with wide wavelength coverage, obtained with FEROS 
at the 2.2m telescope at ESO (La Silla, Chile). 

Effective temperatures T e g are derived spectrophoto- 
metrically from IUE fluxes and Johnson and 2MASS 
magnitudes for all the stars except HR2928 where T e g 
is in agreement with Kilian (1992). Further constraints 
can be derived from the Hei/n ionization equilibrium 
for rSco and HR3055. Figure ^ displays best fits to the 
spectrophotometry. Gravities log g are derived from line 
profile fits to US, H/3, H7 and Ha. Examples for some H 
and He 1/11 lines in r Sco and HR 5285 are given in Fig. [3 

After constraining T c ff and log g we compute small 
grids of synthetic spectra for different carbon abundance 
e(C)= log(iVc/-/VH) + 12 and several values of micro- 
turbulence £. These grids are used for C abundance 
determinations from observation via a ^-minimization 
technique. The free fitting parameters are e(C), vsmi 
and macroturbulence (. The macroturbulent velocities 
in the sample stars remain smaller than twice the sound 
speed. The microturbulent velocity is fixed in the usual 
approach, by demanding e(C) to be independent of the 
equivalent width Wa of the C II line ensemble (see Fig.[5J). 

An excellent match between theory and ob- 



servation is achieved, as shown in Fig. [21 ex- 
emplarily for the hottest (r Sco) and coolest 
(HR5285) star in our sample. The complete 
linelist includes the Cn AA3919.0/20.6, 4267.0/. 2, 
5133.3/37.3/39.2/43.4/45.2/51.1, 5648.1/62.5, 6151.5, 
6461.9, 6578.0/82.9, 6779.9/80.6/83.1/87.2/6791.5 and 
6800.7 A transitions. Fits to Cn A6462A for all the 
stars are shown in Fig. This line together with C II 
A615lA is subject to marked non-LTE effects, turning 
from absorption at spectral type B2 into emission at 
earlier spectral types, a behaviour never reported before. 
The emission results from a non-LTE overpopulation 
of the upper levels of these transitions relative to the 
lower levels, facilitated by close collisional coupling of 
the former to the Cm ground state, which is close to 
LTE in all the sample stars. For HR3055 and HR2928 
Cn AA6462 and 6151 A are not considered in deriving 
the average e(C). See S96 for a discussion of the nature 
of the non-LTE effects of Cn AA4267 and 6578/82 A. 

Non-LTE and LTE abundances for all individual lines 
are displayed as a function of W\ in Fig.0 showing excel- 
lent consistency in non-LTE. A slight degradation of the 
overall consistency is indicated for r Sco (see below) . At- 
mospheric parameters and averaged e(C) are also given. 

4. DISCUSSION & RESULTS 

We have investigated the sensitivity of the transitions 
to modifications of atmospheric and several atomic pa- 
rameters qualitatively. The individual lines show a differ- 
ent behaviour: Cn AA6151 and 6462 A are very sensitive 
to changes of log g, £ and collisional and photoionization 
cross-sections, Cn A4267A is sensitive to photoioniza- 
tion, and all C 11 line strengths react sensitively on T e g 
changes. As an example, exchanging our current photo- 
ionization data with the homogeneous set of N02a would 
result in considerably reduced uncertainties in the mean 
e(C) of t Sco and HR 3055. However, for the same model 
configuration abundances from C 11 A4267 A in the other 
stars are reduced by up to 0.5 dex, approximately repro- 
ducing the LTE results in Fig. while the remaining 
lines behave similarly to the final model. Further analy- 
ses have to be made to quantify these dependencies. 

LTE analyses may produce abundances from the 
prominent Cn A4267A transition in error by ^0.3- 
0.8 dex, by up to ~0.4dex in the case of Cn 
AA6578/82A(note that ~zero corrections may also oc- 
cur), and by ^0.6-0.8dex for the weak Cn AA6151 and 
6462 A lines in the cooler stars. In the hotter stars the 
last two lines turn into emission and cannot be repro- 
duced at all assuming LTE. All other transitions are sub- 
ject to non-LTE corrections on the order of ~0-0.2dex. 

Noteworthy (small) discrepancies to the overall ex- 
cellent agreement arise only for the C 11 AA4267 and 
6578/82 A transitions in rSco. Among the sample stars 
t Sco is the only object showing a considerable (clumped) 
stellar wind and hard X-ray emission (see Howk et 
al. 2000). The problems may be related to these com- 
plications, as wind emission affects the Ha profile and 
the X-ray emission the photoionization rates, which both 
have to be modelled correctly in order to reproduce these 
strongly non-LTE affected C 11 lines accurately. 

Our non-LTE computations reproduce the C 11 
A4267 A theoretical equivalent widths of S96 (his Fig. 1). 
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For Cn AA6578/82A we reproduce the S96 values at 15 
kK; however at 20 kK our W\ are ~10% lower and at 30 
kK up to 50% higher (see S96, Figs. 5 and 6). The dif- 
ference arises because we use non-LTE populations when 
computing the Ha line opacities which define the contin- 
uum against which these lines are measured (S96 assumes 
LTE). Note also that the T e ff-scale employed by S96 for 
his comparison with observation appears to produce val- 
ues of T e ff lower by 700 -3 000 K than our derivations. 

A highly consistent mean e(C) = 8.29±0.03 is derived 
from the sample stars, which provides a tight estimate 
to the present-day C abundance from young stars in the 
solar vicinity. The atmospheric composition appears to 
be unaffected by chemical mixing in the course of stellar 
evolution, i.e. we find no trend of e(C) with evolution- 
ary age. For comparison, adopting Kilian's (1992) results 
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Fig. 1. — Best fits of theoretical energy distributions to measurements by IUE (dotted lines) and Johnson and near-IR 2MASS photometry 
(diamonds). The observed fluxes are dereddened by the values indicated using a standard reddening law and assuming Ry = Ay/E(B — 
V) = 3.1 as typical for the local ISM. They were degraded to the resolution of the ATLAS9 fluxes. The models are normalized to the 
observed V magnitudes and shifted for clarity relative to each other. See Fig. |5]for atmospheric parameters. 
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Fig. 2. — Comparison of synthetic H and Hel/ll lines (smooth lines) with observation of a B0V (r Sco) and a B2V (HR5285) star. 
Atmospheric parameters as summarized in Fig.|S]are employed. 
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Fig. 3. — Examples of best fits (smooth lines) to observed Cn features for a BOV (rSco) and a B2V (HR5285) star. Atmospheric 
parameters and averaged e(C) of each star are given in Fig. 151 
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Fig. 4. — Best fits to the Cn A6462A line in the six sample stars. For clarity, the stronger telluric features in this spectral region are 
suppressed in the observed spectra (dots). 
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Fig. 5. — Non-LTE (filled circles) and LTE (open circles) abundances vs. equivalent width for all the lines that could be measured in 
each spectrum (a linelist is provided in the text). In the upper left part of each row an identification of the star and its carbon abundance 
is given. In the lower left part T e fj, logg, microturbulence, vsini and macroturbulence (where non-zero) are shown. For abundances and 
velocities statistical la- uncertainties are provided and for T c g and logg the errors from our derivation are presented. The grey rectangles 
correspond to lex-uncertainties of e(C). Identification of lines with high sensitivity to non-LTE effects is displayed. Emission lines are 
marked by crosses (C II AA6151 and 6462 A in tSco and A6462 A in HR1861): LTE calculations are not able to reproduce them. 



